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Die Erfindung bezieht sich auf ein 
iteratives Verfahren zur experimentellen Ermittelung der Parameterwerte fiir ein optimales Prozefiergebnis 
innerhaib von vorgegebenen Bereichsgrenzen fur eine Vielzahl von Parametern nicht bekannter Interdepen- 
danz und Ergebnis-Relevanz durch Aufeinanderfolge von Versuchsreihen mit jeweils innerhaib der Bereichs- 
s grenzen abgewandelten Parametervorgaben bzw. auf ein Verfahren zur parameter-abhangigen Ergibnis- 
Kartographierung. 

Dabei sollen insbesondere innerhaib komplexer Systeme mit nicht mehr uberschaubaren oder unbe- 
kannten Abhangigkeiten OptimaParameter mit moglichst geringem Aufwand an praktischen Versuchen 
ermittelt werden. Es konnen ein oder oder mehrere Ziele angestrebt werden. Als Beispiel fur ein 
10 erfindungsgemaB zu losendes Problem ist etwa die Ermittlung der Nahrstoff-Zusammensetzung eines 
Fermentationsmediums zur Bildung eines bestimmten Produktes mit optimaler Raumzeitausbeute zu 
nennen, ggf. unter Beachtung minimaler Herstellungskosten. 

Biotechnologische Prozesse sind exemplarische Beispiele fur schwer modellierbare und experimentell 
relativ gut beherrschbare Systeme. Eine weitere Eigenschaft sind die stark verrauschten Daten der 
75 Experimente. 

Eine grundlegende Rolle bei der Entwicklung biotechnologischer Prozesse spielt nach der Stamment- 
wicklung (Isolierung, Mutation, genetischer Modifikation und Selektion), die Optimierung der Reaktionsbedin- 
gungen. 

Neben pH-Wert und Temperatur ist dabei die qualitative und quantitative Zusammensetzung des 
20 Nahrmediums von entscheidender Bedeutung. 

Die fur jeden Stamm individuelle optimale Nahrmedien-Zusammensetzung beinhaltet ein hohes, wenig 
genutztes Optimierungspotential und ist ein wichtiger okonomischer Faktor eines Fermentationsprozesses. 

Im allgemeinen werden neben der aufwendigen und wenig effizienten Variation nur eines Medienbe- 
standteils ("one-faktor-at-a-time" Methode) Methoden der statistischen Versuchsplanung eingesetzt, um die 
25 Zahl der Versuche zu verringern. In vielen Bereichen, beispielsweise in der ProzeSentwicklung bei 
chemischen Verfahren oder der Qualitatssicherung, wird sie angewandt. 
Problematisch fur ihren Einsatz in der Biotechnologie sind: 
exponentielle Zunahme der Versuche mit der Parameter-Zahl, 
Annahme eines linearen Modells, 
30 Empfindlichkeit gegenuber MeBrauschen. 

Bei funf Parametern und zwei Konzentrations-Stufen benotigt man 2 5 = 32 Experimente, fur zehn 
Parameter schon 2 10 = 1024 Versuche. Will man eine feinere Rasterung der Parametereinstellungen 
vornehmen, z.B. funf Stufen bei funf Parametern, benotigt man 5 s = 3125 Experimente. Der hohen 
Versuchsanzahl kann man durch Verkleinerung der Versuchsplane entgegentreten, was wiederum zu 
35 Vernachlassigungen der Wechselwirkungen zwischen den einzelnen Parametern fuhrt. Die Annahme linea- 
rer Modelle ist gerade bei den nicht linearen biologischen Systemen eine ungunstige Vereinfachung. 

Das biologischen Systemen eigene Problem der schlechten Reproduzierbarkeit der Daten unter nicht 
standardisierten Bedingungen wirft ein weiteres Problem auf. Mit statistischen Methoden kann man das 
Rauschen erfassen und so eine Aussage uber die Qualitat der Daten bekommen. Letztendlich wird man 
40 aber bei diesen Optimierverfahren Messungen mit prozentualen Standardabweichungen oberhalb von 5-10 
% wiederholen mussen. 

Auch andere Optimierverfahren, wie z.B. die Rosenbrock- oder die Simplex-Methode, scheitern in 
komplexen, hochparametrigen Systemen aufgrund der bei ihrer Entwicklung zugrundegelegten Vorstellun- 
gen der zu untersuchenden Systeme (siehe Schwefel H.-P. (1977), Numerische Optimierung von Computer- 
45 Modellen mittels der Evolutionsstrategie, Birkhauser Verlag, Basel und Stuttgart). 

Ziel der Erfindung ist daher eine Versuchsstrategie, die mit einer moglichst geringen Anzahl von 
Versuchen auskommt und trotzdem eine weitgehende Sicherheit bietet, daB Optima "nicht Ubersehen" 
werden. 

Das zu diesem Zweck entwickelte erfindungsgemafie Verfahren der eingangs genannten Art ist dadurch 
so gekennzeichnet, daB die Parametervorgaben jeder Versuchsreihe - ausgehend von den Ergebniswerten 
einer ersten, durch Zufallswerte der Parameter festgelegten Versuchsreihe - durch Anwendung einer dem 
genetischen Evolutions-Algorithmus (nach Holland) folgenden Optimierungsstrategie unter Zugrundelegung 
der Ergebniswerte der jeweils vorangehenden Versuchsreihe festgelegt werden. 

Eine solche Optimierungsstrategie fur die Auswahl von Parameterwerten fur experimentell durchzufuh- 
55 rende Versuchsreihen bei der am Ergebnis einer Versuchsreihe gewichtete, zufallsmodifizierte Zahlenwerte 
Basis fur die Parametervorgaben der nachstfolgenden Versuchsreihe sind, ist ntitzlich bei der experimentel- 
len Suche nach Optimierungen innerhaib von Systemen mit einer sehr groflen Vielzahl von EinfluBgroBen 
(Parametern), deren Abhangigkeit untereinander und Einwirkung auf das zu erzielende Ergebnis weitgehend 
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unbekannt sind oder auch wenn mehrere Zielsetzungen gleichzeitig zu beachten sind, urn mit erheblich 
vermindertem Versuchsaufwand zu Optimierungen zu gelangen. AIs "ProzeBergebnis" wird dabei ein durch 
direkte Einflufinahme auf die praktischen Voraussetzungen veranderbares Ergebnis" verstanden. 

Die erfindungsgemafle Versuchsplanung kann einfach unter Anwendung des Optimieaingsprinzips nach 
5 Holland durchgefuhrt werden, selbstverstandlich sind jedoch auch die zwischenzeitlich in grofler Vielfalt 
entwickelten Abwandlungen desselben einsetzbar, wie sie zusammenfassend von D.E.Goldberg in "Genetic 
Algorithms in Search Optimization and Machine Learning", Addison Wesley Publishing Company, Inc. 
Reading, Massachusetts, 1989 sowie von Lawrence Davis "Genetic Algorithm Tutoriar in Handbook of 
Genetic Algorithms (ed. L.Davis), Van Nostrand Reinhold Verlag New York, 1991 dargestellt sind und 
w insbesondere die Abwandlung der Einzelschritte z.B. in Form von Mehrpunkt-crossover, oder unterschiedli- 
cher Selektionsmechanismen oder auch die Angabe von "Hybriden Algorithmen" betreffen. 

Eine zweckmafiige Realisierung des erfindungsgemafien Verfahrens ergibt sich dadurch, dafi 
la) eine der beliebig wahlbaren Anzahl der Versuche der primaren Versuchsreihe entsprechende Anzahl 
von Bitstrings mit einer der Anzahl der Parameter entsprechenden Anzahl von Segmenten mit jeweils der 
75 gewunschten Rasterung entsprechenden Bitzahl durch Zufallskodierung gebildet wird, 

lb) deren numerisches Ergebnis multipliziert mit dem Quotienten aus Abtastbereich durch maximalen 
decodierten Bitwert zum Parameter-Minimalwert addiert und ggf. aufgerundet je Segment die Realwerte 
der einzelnen Versuchsparameter der primaren Versuchsreihe ergibt, dafi 

II) die Versuche der primaren Versuchsreihe mit den so ermittelten Parametern durchgefuhrt werden 
20 unter Erzielung einer Gruppe von Ergebnissen und dafl 

III) aus den erhaltenen Ergebnissen der Ergebnismittelwert gebildet wird, der einen Wichtungsschlussel 
(aus Einzelergebnis der Versuchsreihe / Ergebnismittelwert der Versuchsreihe) liefert fur die Zufallsaus- 
wahl einer intermediaren Bitstring-Foige, die nacheinander mit- abgestufter Wahrscheinlichkeit einem 
Cross-over (a), der Inversion (b) und der Mutation (c) unterworfen wird, wobei 

25 a) innerhalb der Bitstring-Folge eine Zufallspaarung uber alle Bitstrings vorgenommen wird und bei 

den so gebildeten Bitstringpaaren ein Austausch der einer zufallig gewahlten Position folgenden Bits 
untereinander vorgenommen wird, worauf 

b) hinter einer oder zwischen mehreren zufallen Positionen der mit gegenCiber (a) mit geminderter 
Wahrscheinlichkeit ausgewahlten Bitstrings eine Umkehr der Bit-Reihenfolge unter Beibehaltung ihrer 

30 Parameterzugehdrigkeit und Platznummer stattfindet und schliefilich 

c) mit wiederum verminderter Wahrscheinlichkeit ein Zufallsaustausch einzelner Bits vorgesehen wird, 
dafi die auf diese Weise erhaltene Bitstringfolge der anschliefienden Versuchsreihe gemaS 

lb in Realwerte umgesetzt wird, wobei ggf. der gemaB Wichtung beste Bitstring unverandert in die nachste 
Bitstringfolge eingereiht wird; und dafi nach Durchfuhrung der Versuche entsprechend der so festgelegten 
35 Versuchsvorschriften ggf. wiederkehrend erneut wie unter III verfahren wird. 

Bei der vorstehenden Definition der Erfindung wurden zur Erleichterung des Verstandnisses "Bits" und 
"Bitstrings" fur Strukturelemente und Strukturelementketten eingesetzt, die selbstverstandlich auch durch 
jedes andere Zahlen- oder Buchstabensystem realisiert werden konnen. Auch solche Realisierungen sollen 
hier miterfaflt werden. 

40 Die Genetischen Algorithmen (GA) wurden von J.H. Holland 1975 publiziert ("Adaption in natural and 
artificial systems", The University of Michigan Press Ann Arbor, Michigan). Bei den hier eingesetzten GA 
wurden auch Methoden integriert, die von anderen Autoren stammen. Diese Techniken wurden ursprunglich 
mit dem Ziel erarbettet, eine Simulation der Entwicklung von Systemen mit kompiexen Abhangigkeiten auf 
einem Rechner zu ermoglichen. Beispiele sind: 

45 Funktionenoptimierung; 
Bilderkennung; 
Spieltheorie; 

Simulation der Evolution; 
Parameteridentifikation; 

so Auslegen von Schaltkreisen; 
Optimierung von Gaspipelines. 

Erfindungsgemafl wird nun dieser bislang als Rechenmodell fur EDVisierte Abschatzungen bzw. 
Berechnungen von Entwicklungen, Zusammenhangen oder Verhaitensweisen innerhalb komplizierter Syste- 
me benutzte Algorithmus fur die zweckgerichtete Auswahl von Versuchsreihen zur Erfassung bzw. Katalogi- 

55 sierung der von einer Vielzahl von Parametern abhangenden Ergehnisse und insb. zur Prozefioptimierung 
angegewandt und so eine erhebliche Minderung der Anzahl der durchzufuhrenden Versuche dadurch 
erreicht, dafi ein Vermutungsbereich fur eine Vielzahl von Parametern nicht mehr systematisch in einer 
Vielzahl von Versuchsreihen abgerastert wird, sondern, ausgehend von zufallig gewahlten Punkten innerhalb 
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des Vermutungsbereichs sich abzeichnende Verbesserungen, mit starkerem Gewicht in die nachste 
Versuchsgruppe eingehen, der zusatzlich jedoch wiederum gewisse gewillkurte oder Zufallsanteile beige- 
ordnet werden. 

Der erfindungsgemafie Einsatz von GA ist fur die Optimierung beliebiger Systeme geeignet, fur die eine 
s Abhangigkeit zwischen den variierten Parametern und dem Optimierziel besteht Einige Beispiele sind hier 
angegeben, die aber keinen Anspruch auf Vollstandigkeit erheben. Systeme, bei denen eine Optimierung 
mit GA vorteilhaft sein konnte, sind: 

Zusammensetzung chemischer Reaktionssysteme: 

10 

- Konzentration der Reaktanden 

- Konzentration der Reaktanden und des Katalysators 

- Art des Katalysators, Kombi nation und Katalysatoren 

- Zusammensetzung des Katalysators 

75 - Wahl oder Kombination von Losungsmitteln 

Reaktionsbedingungen chemischer Reaktionssysteme: 

- Temperatur 
20 - pH 

- Druck 

- Reaktionszeit 

- Mischzeit 

- Zeit fur die Zugabe einzelner Reaktanden zum Reaktionssystem 

25 

Biologische Systeme: 

- Zu den Parametern chemischer Systeme kommt noch hinzu: Zusammensetzung des Nahrmediums 

mineralische Medien 
30 komplexe Medien 

Suplementierung komplexer Medien, mit einzelnen definierten Bestandteilen, z.B. Spurenelemen- 
ten oder Vitaminen 

- Begasung 

- Umsetzungen mit Enzymen 

35 

Betriebsweise analytischer Systeme am Beispiel der HPLC: 

- Mischungsverhaltnis der Losungsmittel 

- Art des Gradienten 
40 - Saulentemperatur 

- Druck 

- Flu8 

Optimierung der Rezeptur fur Farben, Aromen, Geruche etc. 

Die Anzahl der Versuche der primaren Versuchsgruppe bzw. -reihe, welche die Zahl der Bitstrings gibt, 
45 wird zweckmaBigerweise zwischen 5 und 160, insbesondere zwischen 5 und 30 gewahlt. Zwei Bitstrings pro 
Untersuchungsparameter sind angemessen. 

Die Anzahl der Bits pro Parameter (Segment), durch welche die Feinheit der Rasterung bestimmt wird, 
kann von Parameter zu Parameter unterschiedlich gewahlt werden und grundsatzlich beliebig hoch sein. 
Eine Bitzahl von 1 bis 10 pro Parameter (1 bis 1024 Unterteiiungen des Vermutungsbereichs) erscheint im 
so allgemeinen als ausreichend. Als "Rasterung" wird dabei die Anzahl der Unterteiiungen (in codierter Form) 
verstanden. 

Die Wahrscheinlichkeiten fur Cross-over, Inversion und Mutation werden insbesondere in den Bereichen 
0,25 - 1; 0 - 1 und 0 - 1 gewahlt. 

Zum Cross-over werden je zwei Bitstrings einander zufallig zugeordnet und ein oder mehrere beliebige 
55 Cross-over-Punkte festgesetzt. 

Die Cross-over-Operation bei 4 Bitstrings mit Vorzugswichtung von II gestaltet sich fur einen Cross- 
over-Punkt z.B. wie folgt: 
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I 0110 11^ 01-100 



II 1100 10^7 Cross-over 11001 

II IIIOOO^S ~~ ^ 11011 

IV 10 1011^ 10000 



Wahrscheinlichkeit fur Cross-over liegt zweckmafiigerweise zwischen 0,25 und 1, vorzugsweise zwi- 
schen 0,6 und 1 . 

Durch Inversion wird die Reihenfolge der einzelnen Bits im Bitstring vertauscht. Die Bedeutung der 
einzelnen Bits andert sich dadurch nicht. Dies erreicht man durch eine Numerierung der einzelnen Bits wie 

z.B.: 



vor Inversion nach Inversion 

Position: 12 | 3456 
A = 1011110 



78 Position: 12|6543|78 

ii a = lobnilii 



(die Inversions-Punkte sind durch Striche dargestellt), 

wodurch die Bedeutung eines Bits unabhangig von seiner Position wird. Inversion wirkt sich erst in den 
nachsten Generationen zusammen mit Cross-over aus. Durch diese Trennung von Position und Bedeutung 
eines Bits sind einige Optimierprobleme besser losbar, z.B. und vereinfacht konnen bestimmte Parameter- 
kombinationen durch Inversion im String raumlich naher positioniert und somit durch Cross-over schlechter 
zerstort werden. 

Inversion ist also die Vertauschung der Reihenfolge der einzelnen Bits unter Beibehaltung ihrer 
Bedeutung, an einem oder zwischen mehreren zufalligen Punkten der Bitstrings. Die Wahrscheinlichkeit fur 
Inversion liegt zwischen 0 und 1, vorzugsweise zwischen 0 und 0,2. 

Mutation variiert den Bitstring-Satz durch punktuellen Austausch eines Bits und verhindert dam it 
vorzeitige Konvergenz. D.h., mit einer bestimmten Wahrscheinlichkeit zwischen 0 und 1 (vorzugsweise 
0,0001 - 0,01) wird im Bitstring eine Null gegen eine Eins bzw. umgekehrt ausgetauscht. 

Die entstandene neue Population ersetzt die alte Generation mit einer Wahrscheinlichkeit von G = 0,3 
bis 1 ,0. Ein Anteil von N * (1 - G) der alten Population kommt unverandert in die neue Population. G liegt 
vorzugsweise zwischen 0,8 und 1 ,0. 

Anschliefiend erfolgt segmentweise die Transformation in die Realwerte der Parameter fur die folgende 
Versuchsgruppe. 

ZweckmaBierweise erfolgt erfindungsgemafi eine Umwandlung in zu maximierende Funktionen: In 
normalen Operations Research Arbeiten wandelt man ein Minimierungs- in ein Maximierungsproblem durch 
die Multiplikation mit (-1) um. Dies ist bei Genetischen Algorithmen nicht anwendbar. Bei Genetischen 
Algorithmen ublich ist folgende Transformation zur Oberfuhrung in ein Maximierungsproblem: 

f(x) = C max - g(x) fur g(x) < C max (1) 
f(x) = 0 andernfalls. 

C m ax ist ein Koeffizient, der beispielsweise das grofite g(x) der Population oder der k letzten Populatio- 
nen sein kann. 

Die Transformation in positive Werte wird durch 

f(x) = u(x) + C min fur u(x) + C min > 0 (2) 
f(x) = 0 andernfalls 

realisiert. Auch hier kann C min frei gewahlt werden, z.B. als schlechtester Wert der aktuellen oder mehrerer 
vorheriger Populationen. k liegt zwischen 0 und 7, wobei 0 konstant die 0 Generation, eins die aktuelle 
Generation und 2 - 6 die entsprechend zuruckliegenden Generationen bedeutet. 7 heiBt. dafi die Umwand- 
lung nicht vorgenommen wird. Vorzugsweise liegt k zwischen 1 und 3. Zweckmafiigerweise ist C min der 
kleinste u(x) der aktuellen oder einer vorangehenden Population. 
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Skalierung: Die Guten werden linear skaliert Mit f als transform ierter und f als ursprunglicher Gute ist 
V = at + b. (3) 

a und b mussen so gewahlt werden, daB die durchschnittlichen Guten von f und f identisch sind, so dafi 
durchschnittliche Mitglieder weiterhin genau einen nachfolgenden Versuch erzeugen. Desweiteren kann 
durch den eingefuhrten Koeffizienten C mult die Anzahl der "Nachkommen" des Besten gesteuert werden, 
wenn gilt f max = C mu i t * f avg . (C mu „ liegt bei kleinen Populationen im Bereich zwischen 1,2 und 10, 
insbesondere zwischen 1 ,2 und 2). 

Selektion: ("survival of the fittest") und anschlieBende Reproduktion dienen zur Vervielfaltigung von 
Individuen (Parametersatzen), die eine uberdurchschnittliche Anpassung (Gute) aufweisen. Die Wahrschein- 
lichkeit, "Nachkommen" fur die nachste "Generation" zu zeugen, ist proportional zum Verhaltnis der Gute 
des einzelnen Individuums zum Durchschnitt aller Guten der aktuellen Population (Bewertungszahl r). 

Die verwendete Selektionsstrategie wurde von James E. Baker, "Reducing bias and inefficiency in the 
selection algorithm" in: Genetic Algorithms and their Applications: Proceedings of the Second International 
Conference on Genetic Algorithms, S. 14-21, Lawrence Erlbaum Associates, Publishers Hillsdale, New 
Jersey 1987, vorgeschlagen: 

1 Ermittlung einer Zufallszahl x; (0 £ x £ 1) 

2 Aufsummieren der Bewertungszahl r jedes Individuums bis: Er>x 

3 Auswahl dieses Bitstrings 

4 x = x + 1 

5 Wenn £ r > x dann nochmalige Auswahl des Individuums, sonst weiter bei 1 

6 Wiederholung, bis die gewunschte Anzahl an Bitstrings erreicht ist. 

Eine Optimierungsstrategie nach Pareto "Cours d'Economie Politique", Rouge, Lausanne (1896) zur 
gleichzeitigen Optimierung mehrerer ZielgroBen wird dadurch realisiert, daB die Population in Lfntergruppen, 
gleicher oder unterschiedlicher GroBe, entsprechend der Anzahl der Ziele aufgeteilt wird und fur jede 
Gruppe ein Selektionsverfahren durchgefuhrt wird, wobei jedoch uber alle Mitglieder der Population 
selektiert wird. 

Die so ermittelte intermediare Population wird nun den genetischen Operatoren unterworfen. 
Es folgen Beispiele fur die Durchfuhrung der Erfindung: 

Beispiel 1: 

Experimentelle Ermittlung der Mediumszusammensetzung (5 Komponenten) fur maximale Zellmassekon- 
zentration beim aeroben Totalabbau von Benzoat mit Pseudomonas put/da 

Das aerobe Bakterium Pseudomonas puttda ist in der Lage, in Wasser gelostes Benzoat zusammen 
mit Sauerstoff zu Zellmasse, Energie, Kohlendioxyd und Wasser umzusetzen. Die vollstandige Umsetzung 
von Benzoat kann jedoch nur erfolgen, wenn alle anderen zum Zellmasseaufbau notwendigen Substanzen 
(Hauptelemente wie Stickstoff, Phosphor, Schwefel..., Spurenelemente wie Mg, Mn, Zn... sowie eventuell 
Wuchsstoffe wie Aminosauren, Vitamine...) in der waBrigen Losung ( , Nahrmedium , ) in ausreichendem Mafie 
vorhanden sind. 

Die optimale Zusammensetzung dieses organismusspezifischen Nahrmediums soil (in Bezug auf die 
umzusetzende Benzoatmenge) fur 5 ausgewahlte Mediumskomponenten (siehe Tab. 1) experimentell 
ermittelt werden. Die anderen notwendigen Mediumskomponenten werden in diesem Beispiel konstant 
gehalten (10 mmol/l Benzoat, 28 mmol/l NH 4 CI, 204 umol/l CaCI 2 , 1.5 mmol/l H 3 P0 4 , 0.534 umol/l C0SO4, 
1.2 umol/l CuSO*. 36.4 umol/l H3BO4, 1.24 umol/l Na 2 Mo0 4 ). 

Die experimentellen Untersuchungen (Versuche) werden in 100 ml Schuttelkolben mit 10 ml Reaktions- 
volumen im Batch-Ansatz durchgefuhrt ( f Optimierung der Start-Mediumszusammensetzung'). Die Reaktions- 
bedingungen werden test vorgegeben: 
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- Animpfmenge (gewaschene Zellen aus Lagerung bei -20 


1 % 


• C):exponentieller Wachstumsphase, Lagerung bei -20 * C): 




- Temperatur : 


37 -C 


- Schuttelfrequenz: 


1 50 rpm 


- Prozeflzeit: 


16.5 h 


- MefigroBe: 


OD G6 o 


(Optische Dichte bei einer Wellenlange von 660 nm, 10 mm Kuvette: als 


MaB fur die Zellmassenkonzentration) 





Die zu variierenden Konzentrationsbereiche fAbtastbereiche') der einzelnen Mediumskomponenten 
werden so gewahlt, dafi die Konzentration jeder Mediumskomponente um etwa 1 Zehnerpotenz verringert 
und etwa 1 Zehnerpotenz gegenuber dem 'Ausgangsmedium' erhoht werden kann (siehe Tab. 1). Als 
Ausgangsmedium stand in diesem Fall ein der Literatur entnommenes Medium zur Verfugung. 



Komponente 


Stamm- 


Medium 


Medium 


Medium 


Einheit 




losung 


Minimal 


Aus gang 


Maximal 




KC1 


100 


0 


1.5 


10 


mmol/1 


MgS0 4 


40.6 


0 


0.609 


4.06 


mmol/1 


FeS04 


719 


0 


10.8 


71.9 


^ntiol/1 


MnS0 4 


1180 


0 


17.7 


118 


pmol/1 


ZnS04 


278 


0 


4.17 


27 .8 


jimol/1 



Tab. 1: Konzentrationen, mogliche Konzentrationsbereiche 



der zu variierenden Mediumskomponenten 

Die Rasterung, also die Anzahl der moglichen Konzentrationsstufen im jeweiligen Abtastbereich (in der 
statistischen Versuchsplanung auch als 'Niveaus' bezeichnet) wird mit 200 + 1 (Nullwert) = 201 vorgegeben, 
um die Optimal konzentration jeder Mediumskomponente moglichst genau ermitteln zu konnen. 

Die technische Realisierung dieser Rasterung erfolgt durch Zugabe einer Stammldsung der jeweiligen 
Mediumskomponente. Da die verwendeten Pipetten nur eine minimale Auflosung von 5 n\ gestatten, ergibt 
sich als maximal dosierbares Volumen je Mediumskomponente 1 000 ul (Minimalvolumen: 0 uJ). Die 
optimalen Volumina der einzusetzenden Stammlosungen ( - Konzentrationen der Mediumskomponenten) 
sollen experimentell ermittelt werden: 

Die Codierung der Volumina der zu variierenden Stammlosungen ('Parameter') erfolgt als bitstrings. Um die 
Rasterung von 201 realisieren zu konnen, sind 8 bit pro Parameter notwendig (2 s = 256). Damit ergibt sich 
eine string-Lange von 40 bit, um eine komplette 'Versuchsvorschriff (Mediumszusammensetzung, Mengen 
der einzusetzenden Stammlosungen) zu codieren. 

Die Versuche werden im Parade lansatz ('Generation') durchgefuhrt. Die Anzahl der parallel durchzufuh- 
renden Versuche wird auf 12 festgesetzt, das bedeutet 2 Parallelversuche je Parameter (5 Stammlosungen) 
und je Ziel (maximale Zellmassenkonzentration gemessen als OD). 

Die Parameterwerte der 0. Generation (Start-Generation) werden bis auf die Startmediums-Zusammen- 
setzung in Binarcode-Darstellung durch Zufall ermittelt. Die Decodierung dieser Binarcode-Darstellung in 
das gewohnte decadische Zahlensystem erfolgt wie ublich: 

Beispiel fur einen Parameter : 

00010011 = 0"2 7 +C2 6 +0*2 5 +1*2* + 0*2 3 + 0"2 2 + T2 1 +1*2° = 19 

Die Umrechnung der decodierten Zahl in den Parameterwert (Menge der einzusetzenden Stammlosung, ul) 
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erfolgt wie angegeben: 
Parameterwert = 

numerisches Ergebnis * Abtastbereich / maximaien 
5 decodierten bit-Wert + Parameter-Minimalwert 

analog obigem Beispiel: 

Parameterwert = 19 * (1 000 - 0) ul / 256 + 0 = 74.2 ul. 
io gerundet 75 ul. 

Die 12 Versuche der 0. Generation werden unter standardisierten Bedingungen parallel durchgefuhrt 
und die ZielgroBe (OD) gemessen (Mengen der eingesetzten Stammlosungen und MeBergebnisse der 
Generation 0 siehe Tab. 2). 

75 Zur Auswertung der in der 0. Generation erhaltenen Ergebnisse wird eine Transformation (Umwandlung 
in eine zu maximierende Funktion, Vermeidung von negativen Ergebnissen) und eine Skaiierung der 
MeBwerte vorgenommen: 

Da es sich in diesem Beispiel bereits urn ein Maximierungsproblem handelt, erfolgt zur Transformation 
lediglich die Addition des niedrigsten Mefiwertes ('C m ^ % = 0.05) zu alien Meflergebnissen (siehe Tab. 2). 
20 Zur Skaiierung der MeBergebnisse werden die Konstanten a und b der Skalierungsgleichung: 

f = a " f + b 

ermittelt (f\ f skalierte, bzw. transformierte 'Gute der Zielfunktion' , hier OD): 

25 

a = favg * (G mu |t * 1 ) / (f max " favg) 
b = favg * 0 - a) 

favg (= 0.242) ist der arithmetische Mittelwert, f max (= 0.506) der Maxim alwert der transformierten f Gute der 
30 Zielfunktion' (OD transformiert). 

Die SteuergroBe C mU | t als MaB fur die Verstarkung des besten Ergebnisses wird in diesem Beispiel auf 
Cmuit = 2 gesetzt. Damit ergibt sich fur die Konstanten der Skalierungsgleichung: 

a = 0.917 b = 0.02 

35 

Mit Hilfe der Skalierungsgleichung werden die transformierten MeBergebnisse skaliert ( siehe Tab. 2). 

Die Berechnung der Bewertungszahl r, als Basis zur Durchfuhrung der 'Selektion* der bisher besten 
Medienzusammensetzungen zur weiteren Bearbeitung mit den f genetischen Operatoren' Crossover und 
Mutation, erfolgt nach: 

40 

r = f/favo" 11/12 



Die Multiplikation mit 11/12 erfolgt, da der bit-string (die Versuchsvorschrift) mit dem besten Ergebnis 
45 unverandert in die nachste Generation (Versuchsreihe) ubernommen werden soil, und daher ein Platz 
weniger fur die Selektion vorhanden ist (Dies ist in diesem Beispiel Versuchsvorschrift Nr. 9, siehe Tab. 2). 



50 
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Nr 


KC1 


MgS0 4 


FeS0 4 


MnS0 4 


ZnS0 4 


OD 


f 




f ' 




r 












Ail 




- 


— 




— 




— 




1 


150 


150 


150 


150 


150 


0 . 066 


0 . 


116 


0 . 


126 


0 . 


477 


2 


125 


800 


785 


955 


100 


0.337 


0. 


387 


0 . 


375 


1 

X • 


42 


3 


695 


105 


445 


925 


530 


0 . 142 


0. 


192 


0 . 


196 


0 . 


742 


4 


845 


510 


265 


5 


315 


0 . 144 


0 . 


194 

X "T 


n 


X -7 O 




7 5 


5 


440 


635 


205 


595 


285 


0.147 


0 . 


1 9,7 

X J / 


o 

\J m 


oni 

•£ V/ X 


n 

\j • 


/ox 


6 


625 


505 


580 


515 


255 


0.291 


o 


341 


o 


J J ,J 


1 

X • 


9 6 1 


7 


600 


895 


60 


525 


960 


0 175 


n 


^ •£ J 


w • 


^ ^ o 


n 

u • 


QC£ 
OJO 


8 


515 


705 


350 


405 


635 


0.242 


0. 


292 


0. 


288 


i. 


091 


9 


955 


610 


835 


35 


380 


0 4 56 

\J m "X -J U 




J v/ O 


n 

\J • 


AQ A 
*± O *t 


l 

X • 


O J J 


10 


95 


665 


270 


830 


355 


0.169 


0. 


219 


0. 


221 


0. 


837 


11 


20 


200 


740 


375 


405 


0.05 


0. 


1 


0. 


112 


0. 


424 


12 


635 


660 


160 


920 


635 


0.08 


0. 


13 


0. 


139 


0. 


527 



Tab. 2: Vorgelegte Mediumskomponenten f gemessene Gute OD, 
transformierte Giite f, skalierte Giite f\ Bewertungszahl r 
der Generation 0. 

30 

Die Selektion erfoigt nach Ermittlung einer Zufallszahl x (in diesem Fall x = 0.63) wie in der Beschrei- 
bung angegeben, siehe Tab. 3): 

Die Bewertungszahl von Versuch 1 ist 0.477. Da die Zufallszahl x grofler als diese Bewertungszahl ist, wird 
35 diese Versuchsvorschrift (bit-string) nicht selektiert. Die Bewertungszahl von Versuch 2 ist 1 .42, womit sich 
als Summe der bisherigen Bewertungszahlen 1.898 ergibt. Jetzt ist die Bewertungszahlsumme groBer als 
die Zufallszahl, so daB Versuchsvorschrift 2 selektiert wird. Die Zufallszahl wird im folgenden Schritt um 1 
auf 1.63 erhoht. Sie ist noch immer kleiner als die Bewertungszahlsumme, so dafl Versuchsvorschrift 2 
erneut selektiert wird. Eine erneute Erhohung der Zufallszahl x um 1 auf 2.63 fuhrt dazu , daB x jetzt groBer 
40 ist als die Bewertungszahlsumme, so daB das Verfahren mit Versuch 3 weitergefuhrt wird u.s.w. (siehe 
Tabelle 3). 
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Nr. 


r 


Summe ( r ) 


x, x+1 


Selektion 
x<Summe(r ) ? 


selektiert 
Nr. 


1 


0.477 


0.477 


A f O 

0 . 63 


nein 




2 


1.42 


1 a ft «j 

1.897 


0 . 63 


ja 


2 






1 . 897 


1 . 63 


ja 


2 






1 O Q *7 


Z . OJ 


nein 






n 749 


0 O Q 


Z • OJ 


j a 


3 






z . b j y 


j ♦ b j 


nein 




4 




3.389 


3.63 


nein 




5 


0.761 


4.150 


3 . 63 


ja 


5 






4.150 


4 . 63 


nein 




6 


1.261 


5.411 


4 . 63 


ja 


6 






5.411 


5.63 


nein 




/ 


U .ODD 


o.2b / 


C /TO 
O • OJ 


j a 


7 






O • Z O / 




nein 




Q 
O 






b . b J 


ja 


8 






7-358 


7 . 63 


nein 




9 


1.833 


9.191 


7.63 


ja 


9 






9.191 


8.63 


ja 


9 






9.191 


9.63 


nein 




10 


0.837 


10 . 028 


9 . 63 


ja 


10 






10 .028 


10.63 


nein 




11 


0.424 


10.452 


10.63 


nein 




12 


0.527 


10.979 


10.63 


ja 


12 



Tab. 3: Selektion, Vorgehensweise bei der Selektion 



Nach Abschlufl der Selektion zeigt sich, dafl in diesem Fall die Versuchsvorschriften 2 und 9 doppelt 
selektiert werden, wahrend die Versuchsvorschriften 1, 4 und 11 nicht ubernommen werden. 

Die selektierten bit-strings werden danach jeweiis paarweise einem crossover unterzogen. Die Auswahl 
der Paare erfolgt zufallig, die crossover-Wahrscheinlichkeit ist zu 95 % festgelegt worden. Die Durchfuhrung 
des crossover ist z.B. fur 2 beliebige bit-strings im folgenden dargestellt, wobei der zufallig gewahlte 
crossover-Punkt durch * markiert ist: 

11110001 10101*100 10110101 01111100 10101111 
00011101 01010*111 00101110 10101001 01010111 

neue bit-strings: 
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11110001 10101111 00101110 10101001 Q.1010111 
00011101 01010100 10110101 01111100 10101111 

Nach der crossover Operation werden an den so neu generierten bit-strings Punktmutation durchgefuhrt, in 
dem eine 1 gegen eine 0 Oder umgekehrt ausgetauscht wird. Die Mutation erfolgt zufallig mit einer 
Mutationswahrscheiniichkeit von 0.1 %. Am Beispiel des oben neu generierten bit-strings ist eine Mutation 
dargestellt, markiert durch eine eckige Klammer: 



00011101 01010100 10110101 01[1]11100 10101111 

neuer bit-string: 



00011101 01010100 10110101 01[0]11100 10101111 

Die so neu erzeugten bit-strings werden anschlieBend decodiert und die Versuche wie oben beschrieben 
durchgefuhrt (1. Generation). Dieser Zyklus wird in diesem Beispiel insgesamt 7 mai durchlaufen (Genera- 
tion 0-6). 

Die Einstellungen des Genetischen Algorithmus fur alle Generationen sind rm folgenden noch einmal 
zusammengefafit: 



SteuergroGe zur Skalierung: 




= 2 


crossover Wahrscheiniichkeit: 




= 95 % 


Mutations Wahrscheiniichkeit: 




= 0.1 % 


Inversions Wahrscheiniichkeit: 




= 0.0 % 


Wahrscheiniichkeit fur neue Generation: 




= 100 % 



Die Crossover-Wahrscheinlichkeit wird wesentlich hoher gewichtet, als die Mutations Wahrscheiniich- 
keit, weil durch Kombination der Parameter schneller das Optimum gefunden wird. Inversion verbessert in 
diesem Beispiel die Suche nicht und wird daher zu 0 gesetzt Die Wahrscheiniichkeit fur eine neue 
Generation (Austausch der bit-string Folge) wird zu 100 % gesetzt, da so der Abtastbereich maximal 
abgesucht wird. 

Ergebnis: Durch Einsatz des Genetischen Algorithmus zur experimentellen Ermittlung der Startkonzen- 
tration von 5 ausgewahlten Mediumskomponenten zur Maximierung der Zellmassenkonzentration beim 
aeroben Totalabbau von Benzoat mit Pseudomonas putida konnte mit insgesamt 84 Versuchen die 
Zellmassenkonzentration urn den Faktor 9.4 gesteigert werden (OD s , ar tmedium 0.045, OD 0 p« ima imedium 0.424). 

Die Zusammensetzung von Start- und Optimalmedium ist in Tabelle 4 dargestellt: 
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Ausgangsmedium Optimalmedium 



Einheit 



KC1 



1.5 



8.3 



mmol/1 



MnS0 4 
ZnS04 



MgS0 4 
FeSQ 4 



4.17 



0.609 



17.7 



10.8 



4 . 14 



0.278 



3.49 



60,1 



mmol/1 
jjmol/1 
pmol/1 
/jmol/1 



Tab. 4: Konzentrationen Ausgangs- und Optimalmedium 

Der Vergleich mit herkommlich zur Mediumsoptimierung eingesetzten Techniken der statistischen 
Versuchspfanung zeigt die Uberlegenheit des erfindungsgemaflen Verfahrens zur experimentellen Ermitt- 
lung der Parameter fur optimale Ergebnisse: 

S. Bloor und R.R. England (Enzyme Microb. Technoi. 1_3 (1991) 76-81) setzten zum Beispiel 'factorial 
experimental design* ein, um ebenfalls die Konzentratiori von 5 Mediumskomponenten zu optimieren. Sie 
fuhrten 16 Versuche durch, um lediglich 2 Konzentrationsstufen (Niveaus) jeder Mediumskomponente zu 
untersuchen (+/- 50 %-Wert des Ausgangsmediums). Damit waren bei 5 Parametern 8 Versuche je 
Konzentrationsstufe erforderlich. Mit dem hier entwickelten Verfahren mit weitaus hoherer Auflosung (201 
Konzentrationsstufen oder Niveaus) sind dagegen nur 0,4 Versuche je Konzentrationsstufe erforderlich. 

Die erforderliche Anzahl der durchzufuhrenden experimentellen Versuche kann also, wie dieses Beispiel 
zeigt, gegenuber Techniken der bekannten statistischen Versuchsplanung um den Faktor 20 reduziert 
werden. 



Experimentelle Ermittlung der Mediumszusammensetzung (12 Komponenten) fur maximale Zellmassekon- 
zentration beim aeroben Totalabbau von Benzoat mit Pseudomonas putida 

Ziel einer 2. Optimierung war es, neben der Maximierung der Zellmassenkonzentration durch Optimie- 
rung der Startkonzentration von 12 Nahrsalzen und Spurenelementen gleichzeitig al!e 12 eingesetzten 
Mediumskomponenten in einer Versuchsabfolge zu minimieren. Es sollten also gleichzeitig 13 Optimierziele 
verfolgt werden, daher wurde nach der von Pareto konzipierten Vektoroptimierung vorgegangen. 

Unter denselben experimentellen Bedingungen wie zuvor wird eine Optimierung durchgefuhrt, lediglich 
die Rasterung wird von 201 auf 101 Konzentrationsniveaus je Mediumskomponente verringert. 

Damit sind noch 7 bit pro Parameter erforderlich (2 7 = 128). Die string-Lange ergibt sich bei 12 
Parametern zu 84 bit. 

Die Konzentration der eingesetzten Stammlosungen, sowie der abgedeckte Konzentrationsbereich sind 
in Tabelle 5 zusammengefafit. 

Die Anzahl der parallel durchzufuhrenden Versuche je Generation wird auf 28 festgesetzt, das bedeutet 
2 Versuche je zu minimierende Mediumskomponente und 4 Versuche fur das Ziel Zellmassekonzentration, 
das somit ein etwas starkeres Gewicht erhalt. 



Beispiel 2: 
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Komponente Stamm- Medium Medium Medium Einheit 
losung Minimal Ausgang Maximal 



TslH a P 1 
Via. A \~ X 


7 7 A 


H3PO4 


200 

^ U u 


XCC 1 
J\l~ X 






81 1 




1440 


MnS0 4 


1180 


Z11SO4 


278 


CUSO4 


160 


C0SO4 


71.1 


CaC12 


27 .2 


H3BO4 


4.85 


Na2Mo0 4 


165 



U 


n n 0 q 


U 


1.3 


U 


1 C 
1.3 


n 
U 


u. buy 


0 


10.8 


0 


17.7 


0 


4.17 


0 


1.2 


0 


0.534 


0 


0.204 


0 


0.0364 


0 


1-24 



n 

\j . 


1 fl 7 
x 0 / 


iilVJ X / X 


1 0 




UUIIU X / X 


i n 
X u 




mmoi/ x 


A 

4 . 


Uo 


mmox/ 1 


7 1 


. 9 


pmol/1 


59 


.2 


/imol/1 


13 


.9 


jimol/1 


8. 


01 


pmol/1 


3 . 


56 


pmol/1 


1. 


36 


mmol/1 


0. 


243 


mmol/1 


8. 


27 


/imol/1 



Tab. 5: Konzentrationen, mogliche Konzentrationsbereiche 
der zu variierenden Mediums komponenten 

Zur Auswertung der Experimente analog zum Beispiel 1 mufi fur die zu minimierenden Parameter eine 
Umwandlung in zu maximierende Werte vorgenommen werden, das heifit die Parameterwerte (g) jeder 
einzelnen Mediumskomponente werden von ihrem jeweiligen maximalen Wert C max substrahiert: 

f = Cmax - g fur g < C max ; 
f = 0 fur g > C max . 

Die Transformation und Skalierung erfolgt wie in Beispiel 1 gezeigt. Bei der Ermittlung der Bewertungs- 
zahl r mufi fur jedes Ziel eine Bewertungszahl r ermittelt werden (13 Bewertungszahlen). Bei der Selektion 
mufi dann entsprechend die bit string Folge in eine Gruppe mit 4 bit strings (Bewertungszahl OD) und 12 
Gruppen mit 2 bit strings (Bewertungszahlen je Mediumskomponente) unterteilt werden. 

Diese Optimierung wurde analog dem Beispiel 1 bis zur 7. Generation durchgefuhrt, der Zyklus also 
insgesamt 8 mat durchlaufen (Generation 0-7). 

Ergebnts: Durch Einsatz des Genetischen Algorithmus zur experimentellen Ermittlung der Startkonzen- 
tration von 12 ausgewahlten Mediumskomponenten zur Maximierung der Zellmassenkonzentration und 
gleichzeitiger Minimierung der eingesetzten Mediumskomponenten beim aeroben Totalabbau von Benzoat 
mit Pseudomonas putida konnte mit insgesamt 224 Versuchen die Zellmassenkonzentration urn den Faktor 
28 gesteigert werden (OD startmedjum 0.045, OD 0ptimalmedium 1.272). 

Die Zusammensetzung von Start- und Optimalmedium ist in Tabelle 6 dargestellt. 
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Komponente 


Ausgangsmedium 


Optima lmedium 


Einheit 


(Benzoat 


10 


10 


mmol/1 ) 


NH4CI 


*5 O 

2 y 


22 • 4 


mmol/1 


H3PO4 


1 . b 


7.9 


mmol/1 


KC1 


1-5 


O.J 


mmol/1 


MgSO^ 


A c r\ c\ 
0.609 


3.49 


mmol/1 


r eb(J4 


1 A O 


a a i 


fimol/ 1 


MII0U4 


1 "7 "7 
1 / . / 


A 1/1 


£/mpl/ 1 




4.1/ 


a no 
U . 2 / 0 


jjmoL/ 1 


CUSO4 


1 . 2 


0 . 881 


jimol/1 


C0SO4 


0.534 


1.32 


jjmol/1 


CaCl 2 


204 


639 


pmol/1 


H3BO4 


36.4 


72.8 


jimol/1 


Na2Mo04 


1.24 


1.16 


jimol/1 



25 

Tab. 6: Konzentrationen Ausgangs- und Optima lmedium 



Auch hier zeigt ein Vergieich von herkommlich zur Mediumsoptimierung eingesetzten Techniken der 
30 statistischen Versuchsplanung die Uberlegenheit des entwickelten Verfahrens zur experimentellen Ermitt- 
lung der Parameter fur optimale Ergebnisse: 

Roseiro et al. (Process Biochemistry 27 (1992) 167-175) setzten zum Beispiel 'factorial experimental design* 
nach Plackett and Burman ein, um die Konzentration von 1 1 Mediumskomponenten zu optimieren. Sie 
fuhrten 28 Versuche durch, um lediglich 2 Konzentrationsstufen (Niveaus) jeder Mediumskomponente zu 
35 untersuchen (+/- 20 %-Wert des Ausgangsmediums). Damit waren bei 11 Parametern 14 Versuche je 
Konzentrationsstufe erforderlich. 

Mit dem hier entwickelten Verfahren mit weitaus hbherer Auflosung (101 Konzentrationsstufen oder 
Niveaus) sind dagegen nur 2 Versuche je Konzentrationsstufe erforderlich. Dabei wird jedoch nicht nur 1 
Ziel wie bei Roseiro et al. verfolgt, sondern gleichzeitig 13 Ziele. 
40 Weitere aktuelle Beispiele fur die bisher eingesetzte Technik 'statistische Versuchsplanung' zur Opti- 
mierung der Startmediums-Zusammensetzung: 
R. G. Silveira et al. : J. of Ferment, and Bioeng. 72, 1 (1991), 20-25. 

R. Banerjee, B.C. Bhattacharyya: Biotechnol. and Bioeng. 41 (1993), 67-71. 

P. V. Rao et al.: Process Biochemistry 28 (1993), 391-395. 

45 

Patentanspruche 



1. Iteratives Verfahren zur experimentellen Ermittelung der Parameterwerte fur ein optimales Prozefiergeb- 
nis innerhalb von vorgegebenen Bereichsgrenzen fur eine Vielzahl von Parametern nicht bekannter 
50 Interdependanz und Ergebnis-Relevanz durch Aufeinanderfolge von Versuchsreihen mit jeweils inner- 

halb der Bereichsgrenzen abgewandelten Parametervorgaben, 
dadurch gekennzeichnet, 

dafl die Parametervorgaben, jeder Versuchsreihe - ausgehend von den Ergebniswerten einer ersten, 
durch Zufallswerte der Parameter festgelegten Versuchsreihe - durch Anwendung einer dem geneti- 
55 schen Evolutions-Algorithmus (nach Holland) folgenden Optimierungsstrategie unter Zugrundelegung 

der Ergebniswerte der jeweils vorangehenden Versuchsreihe festgelert werden. 
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2. Verfahren nach Anspruch 1, 
dadurch gekennzeichnet, 
dafi 

la) eine der beliebig wahlbaren Anzahl der Versuche der primaren Versuchs reihe entsprechende 
Anzahl von Bitstrings mit einer der Anzahl der Parameter entsprechenden Anzahl von Segmenten 
mit jeweiis der gewunschten Rasterung entsprechenden Bitzahl durch Zufallskodierung gebildet 
wird, 

lb) deren numerisches Ergebnis multipliziert mit dem Quotienten aus Abtastbereich durch maxima- 
len decodierten Bitwert zum Pararneter-Minimalwert addiert und ggf. aufgerundet je Segment die 
Realwerte der einzelnen Versuchsparameter der primaren Versuchsreihe ergibt, daB 

II) die Versuche der primaren Versuchsreihe mit den so ermitteiten Parametern durchgefuhrt werden 
unter Erzielung einer Gruppe von Ergebnissen und daB 

III) aus den erhaltenen Ergebnissen der Ergebnismittelwert gebildet wird, der einen Wichtungs- 
schlussel (aus Einzelergebnis der Versuchsreihe / Ergebnismittelwert der Versuchsreihe) liefert fur 
die Zufallsauswahl einer intermediaren Bitstring-Folge, die nacheinander mit abgestufter Wahr- 
scheinlichkeit einem Cross-over (a), der Inversion (b) und der Mutation (c) unterworfen wird, wobei 

a) innerhalb der Bitstring-Folge eine Zufallspaarung uber alle Bitstrings vorgenommen wird und 
bei den so gebildeten Bitstringpaaren ein Austausch der einer zufallig gewahlten Position 
folgenden Bits untereinander vorgenommen wird, worauf 

b) hinter einer Oder zwischen mehreren zufallen Positionen der mit gegenuber (a) mit geminderter 
Wahrscheinlichkeit ausgewahlten Bitstrings eine Umkehr der Bit-Reihenfolge unter Beibehaltung 
ihrer Parameterzugehorigkeit und Platznummer stattfindet und schhiieBlich 

c) mit wiederum verminderter Wahrscheinlichkeit ein Zufallsaustausch einzelner Bits vorgesehen 
wird, 

dafi die auf diese Weise erhaltene Bitstringfolge der anschliefienden Versuchsreihe gemaB 
lb in Realwerte umgesetzt wird, wobei ggf. der gemaB Wichtung beste Bitstring unverandert in die 
nachste Bitstringfolge eingereiht wird; und daB nach Durchfuhrung der Versuche entsprechend der so 
festgelegten Versuchsvorschriften ggf. wiederkehrend erneut wie unter III verfahren wird. 

3. Verfahren nach Anspruch 2, 
dadurch gekennzeichnet, 

daB die primare Versuchsreihe 5 bis 30 Versuche umfaflt. 

4. Verfahren nach Anspruch 2, 
dadurch gekennzeichnet, 

daB pro Parameter zwei Bitstrings vorgesehen werden. 

5. Verfahren nach Anspruch 2, 
dadurch gekennzeichnet, 

daB die Zahl der Bitstrings je Versuchsreihe innerhalb der Aufeinanderfolge von Versuchsreihen 
konstant bleibt. 

6. Verfahren nach einem der vorangehenden Anspruche, 
dadurch gekennzeichnet, 

daB die Wahrscheinlichkeiten fur Cross-over (a); Inversion (b) und Mutation (c) zwischen 0,25 und 1 ; 0 
und 1; sowie 0 und 1 gewahlt werden. 

7. Verfahren nach Anspruch 2, 
dadurch gekennzeichnet, 

daB in Stufe III die erhaltenen Ergebnisse im Falle eines Minimierungsproblems in ein Maximie rungs- 
problem nach 

f(x) = C max - g(x) fur g(x) < C max 
f(x) = 0 andernfalls, 

umgewandelt werden, 

und negative Gutewerte in positive Gutewerte nach 
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f(x) = u(x) + C mjn fur u(x) + C min > 0 
f(x) = 0 andernfalls 

transformiert werden, 

so skaliert werden, daB die Ergebnisse so gestaucht oder gestreckt werden, da8 ein durchschnittliches 
Ergebnis seinen Bewertungsfaktor beibehalt und daB der Bewertungsfaktor einen maximalen Wert nicht 
uberschreitet, mit 

f = af + b 

f max ~ C mu it*favg • 

die Wahrscheinlichkeit fur Cross-over zwischen 0,6 und 1 , vorzugsweise 0,9 und 1 liegt, die 
Wahrscheinlichkeit fur Inversion vorzugsweise 0,0 betragt, die Wahrscheinlichkeit fur Mutation zwischen 
0,0001 und 0,1, vorzugsweise zwischen 0,001 und 0,05 liegt, 

zusatzlich zum Selektionsverfahren der Beste eines Zieles in die nachste Generation ubernommen 
wird, und 

die Population in Untergruppen, gleicher oder unterschiedlicher GroBe, entsprechend der Anzahi der 
Ziele aufgeteilt wird. 

Verfahren nach Anspruch 7, 
dadurch gekennzeichnet, 

daB C max das groBte g(x) der aktuellen Population oder einer vorhergehenden Population ist, 
C m i„ der kleinste u(x) der aktuellen Population oder einer vorhergehenden Population ist, und 
C mu it zwischen 1 ,2 und 1 0 liegt. 

Verfahren nach Anspruch 7, 
dadurch gekennzeichnet, 

daB bei der Ermittlung der intermediaren Bitstring-Folge fur jedes Ziel, fur jeden Bitstring ein 
Bewertungsfaktor ermittelt und fur jedes Ziel ein Selektionsverfahren folgendermaBen stattfindet, 

1. Ermittlung einer Zufallszahl x; (0 £ x ^ 1) 

2. Aufsummieren der Bewertungszahlen r jedes Individuurns bis Summe (r) > x; 

3. Auswahl dieses Bitstrings; 

4. x = x + 1 ; 

5. Wenn Summe (r) > x nochmalige Auswahl dieses Bit-Strings, zuruck zu 4, sonst zu 1 ; 

6. Wiederholung, bis die gewunschte Anzahi an neuen Bitstrings erreicht ist, 
wonach anschlteBend Cross-over, Inversion und Mutation erfolgen. 
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(57) Der experimentelle Aufwand zur Ermittlung der 
Daten fur die Optimalbedingungen eines von einer Viel- 
zahl von EinfluGgrdGen abhangenden Prozesses wird 
durch eine programmgesteuerte Versuchsstrategi e 
stark vermindert, indem ausgehend von Zufallswerten 
fur die einzelnen Parameter und den damit erzielten 
Ergebnissen (erste Versuchsgruppe) eine dem Evoluti- 
ons-Algorithmus folgende Auswahl der zweiten und fol- 
genden Versuchsgruppen von Parameterwerten mit 
Wichtung am Ergebnismittelwert getroffen wird. mit 
zusatzlicher Zufallsabwandlung, wodurch sichergestelrt 
wird, daG vorhandene Optima nicht "ubersehen" wer- 
den. Bei der Realisierung erfolgt eine Codierung der 
Parameterwerte in Form von Bitstrings. deren Anzahl 
der Versuchszahl einer Versuchsgruppe errtspricht und 
die in eine der Zahl der Parameter entsprechende 
Anzahl von Segmenten mit einer der gewOnschten 
Rasterung entsprechenden Brtzahl unterteirt sind. Nach 
jeder Versuchsgruppe erfolgt zur Festlegung der Daten 
der nachsten eine am Ergebnis gewichtete Zufallsaus- 
wahi unter abgestufter Anwendung von Cross-over-, 
Inversions- und Mutationsauswahlprinzipien. 
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